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ABSTRACT 



We present theoretical analysis on old white dwarf stars with some hydrogen, that 
possess a mass of surface hydrogen from 1 x lO~ n M0 to 1 x 1O~ 7 M0. The evolution of 
such objects is complicated by convective mixing from surface convection zone to the 
underlying helium layer. In this paper, we provide first self-consistent, quantitative in- 
vestigation on the subject of convective mixing. Numerical cooling curves and chemical 
evolution curves are obtained as a function of white dwarf mass and hy drogen content 



Such results will be applied to the investigation of the non-DA gap of I Bergeron et al 
(|1997l ) in a later paper. 



Subject headings: (stars:) white dwarfs, stars: atmospheres, stars: evolution, convec- 
tion, diffusion 



Introduction 



White dwarf stars are the burnt-out relics of lower mass main sequence stars and are by far 
the most abundant stellar remnants in the uni yerse. It is bel i eved that up to 97% of the stars 
in our galaxy will end up as white dwarfs (e.g. iFontaine et al.l (120011 ) ). As such, studies of white 
dwarf populations provide invaluable information of the history of star formation and evolution in 
our galaxy. The modeling of white dwarf populations rests on the physics of white dwarf cooling, 
requiring a detailed treatment of matter under a wide range of pressure and temperature. Most of 
the heat capacity resides in the dense carbon/oxygen core. The rate of energy loss, on the other 
hand, is controlled by the thin non-degenerate envelope composed of hydrogen or helium. 

The gravitational acceleration in white dwarf atmosphere is very strong, at the order of 
10 8 m/s 2 . The gravitational field induces a non- vanishing electric field that separates chemical 
species of different charge to mass ratio. In the absence of convective instability, the outermost 
part of a white dwarf atmosphere has a simpli fied structure consisting of an almost pure hydrogen 



layer on top of an almost pure helium layer (jDAntona &: Mazzitellil Il990l : IFontaine et al 
Hansen! 12003 ). 
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The rate of white dwarf cooling is significantly affected by the surface composition (e.g. lHansen 

(|1999l )). It is most straight-forward to assume that the composition of the white dwarf photo- 



sphere is determined by what remains after the mass loss on the Red Giant and Asymptotic Giant 
Branches, i.e., DB and DC stars are different from DA by lacking hydrogen. However, observations 
suggest that the story may be more complicated. 

It is observed that the spectral composition and effective temperature are correlated. Two 
intervals in T e ff space have been found to be almost devoid of non-DA stars. Various mecha- 
nisms of spectral evo lution have been explored and discussed qualitatively in the literature (e.g. 
Bergeron et al.l (|200ll )). The first gap resides at 40000-ftT to 35000-ftT and is interpreted to be the 
result of diffusion completion, follo wed by convective mixing from the underlying helium convection 
zone ([Fontaine Wesemaelill987l ). The explanation o f the other gap (a.k.a. the "non-DA gap"), 



located at 5000K to 6000-ftT, remains rather unclear ([Bergeron et al.l 119971 . 1200ll ). Each of these 
gaps indicates that the photosphere composition of an individual star could be variable during the 
course of white dwarf cooling. The mainstream opinion is that, since a white dwarf can do nothing 
but cool, the two gaps suggests that a white dwarf could change its surface chemical composition 
from helium to hydrogen to helium to hydrogen to helium as it cools, due to the interplay of ac- 
cretion of interstellar hydrogen and convective mixing. However, such statement lacks quantitative 
verification. 

In this series of paper, our goal is to 

1. Quantitatively investigate the interplay between spectral evolution and white dwarf cooling. 

2. Explain the existence of non-DA gap and explore its cosmological implications. 



In the present paper, we hold the opinion that convective mixing (specifically, convective mixing 
from surface hydrogen layer to the underlying helium layer) is responsible for spectral evolution. 
Our statement will be quantitatively validated. We restrict ourselves to consider only the simplest 
scenario of convective mixing without any accretion from the interstellar medium (i.e., a model 
where total hydrogen mass, ma, is fixed.) It will be shown that the cooling curves produced from 
such scenario are already consistent with observations. In principle, accretion can add a further 
dimension to the problem, but is unnecessary to capture the basic features of observation. 

The paper is outlined as follows: In § [2] we briefly review past investigations on convective 
mixing, which is rather controversial. In § [3] we study white dwarf envelopes in which convective 
mixing is absent. We then introduce our self-consistent treatment of convective mixing in §2] and 
revisit the convective mixing scenario in a quantitati ve fashion in jj [5l In section §[6] we incorporate 
our results with the white dwarf evolution code of lHansen! (|1999l ) to calculate numerical cooling 
curves and chemical evolution curves. Finally, we discuss and conclude our results in §[7] and §[HJ 
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2. Brief Review on Convective Mixing 



Here, we specifically refer convective mixing to that from the surface hydrogen conve ction 
zone to the underlying heli um layer. Such scenario was first proposed by IShipmanl (|1972l ) and 
was further investigated by iBaglin fc Vauclairl (|1973l ) and iKoesterl (|1976l ). It can be shown that, 
at most observable T e ff, the surface hydrogen convection zone grows as T e ff decreases, (see e.g., 
fig- E]) Thus, in the proposed scenario, a white dwarf with an intermediate amount of hydrogen (i.e., 
optically thick while less than O(10 -7 M©)) will appear as a DA white dwarf until the base of the 
surface convection zone reaches the underlying helium layer. When the base of the convection zone 
reaches the underlying helium layer, Shipman claimed that the convective motion will dredge-up the 
underlying helium to the surface, resulting in a DB (helium lines) or DC (no lines) white dwarf. 
Shipman also made an order of magnitude estimation (based on the difference in the radiative 
opacity of H and He gas) and claimed that a 13000i"C DA white dwarf would evolve into a 19000-ftT 
DB white dwarf. Such estimation is apparently crude since the existence of convection zone is 
completely ignored. 



Later, IBaglin &: Vauclairl (119731 ) re-investigated the problem by considering models of convec- 
tive envelope. By assuming that the mixing operates at a timescale much longer than hydrody- 
namical timescale and comparing temperature of different envelope models at q = 1 — Jf r = 10 -6 , 



they concluded that the T e ff of the resultant helium white dwarf cannot exceed 18000-fT, i.e., the 
hottest DB stars cannot be evolved from DA stars. 



Koesterl (119761 ) considered convective mixing in his own white dwarf models and showed no 
increase in T e ff upon convective mixing. He claimed that there should not be an increase in T e ff 
because all envelope solutions converge to the s o-called "radiat ive zero solution" independent of the 
boundary conditions. Since the publication of IKoesterl (|1976l ). the opinion of the majority seems 
to be that upon convective mixing, a DA white dwarf will turn into a DC white dwarf at the same 

T eff- 

Despite their disagreement in the evolution of T e ff, all of the above numerical works assume 
that convective mixing would instantaneously evolve the surface from pure hydrogen composition 
to nearly pure helium composition. Such unjustified imposition is probably encouraged by the fact 
that convection zone in pure helium model is generally, orders of magnitude more massive than 
that of hydrogen model. However, as we will show later in this paper, such simplification is not 
always appropriate. 



Vauclair &: Reissd (119771 ) took a even closer look at the convective mixing picture by investi- 



gating the effects of /x-barrier. Strictly speaking, the white dwarf envelope is not exactly layered, 
since classical diffusion would counteract gravi tation-induced diffusion at the H/He layer inter- 
face (jD'Antona Mazzitellilll990l : iHansenl 12004 ) . Such counteraction would result in a /i- gradient, 
which stabilizes the convective envelope. However, It is shown by the authors that /i-barrier cannot 
stop convective mixing as long as 



Mwd 



> 10 



-14 
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On the observation front, iBergeron et al.l (|1990l ) analyzed the high Balmer lines of colder DA 
white dwarfs and found the lines has a tendency to broaden as T e ff decreases. They interpreted 
this result as an evidence of convective mixing. Helium gas is spectroscopically invisible at the T e f f 
of their concern (< 11500i^), hence convective mixing could preserve the original spectral type of 
a white dwarf (DA) while reducing the surface opacity, resulting in a higher photospheric pressure 
and broader spectral lines. 

Thus, the outcome of convective mixing is rather controversial. Shall T e ff increase after the 
mixing? Shall spectral type change? We will address these questions in the following sections. 



3. Structure of White Dwarf Envelope without Convective mixing 

Before we investigate the effects of convective mixing, let us briefly study the structure of white 
dwarf envelope with fixed surface composition (i.e., where convective mixing is absent). A careful 
study of these relatively simple cases yields insights on the topic we will later address. We consider 
two examples: The first being a pure helium white dwarf envelope with My/D = 0.6M & . The other 
is a model with identical Mwd but has a layer of pure hydrogen (run = 10~ 5 M & ) on top of the 
helium layer. As we described in § [TJ this is the equilibrium configuration of white dwarf envelope 
in the absence of convective mixing. 

The method we c onstr uct our envelope model is standard and similar to that described in 
Fontaine Sz van Horn (jl976l ). Ho wever, we men tion a few difference: First, at low T e ff, we use 
realistic radiative transfer code of iHansenl (|1999l ) for surface boundary conditions. This is a major 
i mprovement because envelope structure is very sensitive to the boundary condition at lower T e ff 
([Fontaine et al.Nl974h. Secondly, we implemente d the up-to-date EOS ( equation of state) table 



of 



Saumonetal 



(]1995l ) and the opacity table of iRogers fc Iglesiad (|1992l ). Incidentally, the EOS 
of Saumon et al. as well as our surface boundary condition subroutine has the ability to deal 
with H-He mixture. This feature is not exploited in the current section, however, it is of crucial 
importance as we proceed into section [H 

We calculate the structure of the envelope on a T e ff grid ranging from \ogT e ff = 4.5 to 
log T e ff = 3.35 (all physical quantities in this work are expressed in cgs units unless mentioned 
otherwise). A few representative models are shown in figure [TJ Our results recover the well- 
known fact that white dwarf envelopes generally feature a surface convection zone on top of radia- 
tive/conductive zone which connects to the degenerate core. The relation between surface convec- 
tion zone mass (m cz ) and T e ff f° r hydrogen models are calculated and shown in figure [2} Following 
the cooling sequence, convection zone first deepens with the decrement of T e ff. The curve then 
flattens and gradually inflected at lower temperatures (T e ff < 4000-fT). Note that the maximum 
mass of convection zone is ~ 1.5 x 1CP 6 M©, lower than our assumed hydrogen mass 1(T 5 M . Thus, 
convective mixing needs not to be considered in this example. 



We investigate the relation between T e fj and core temperature T c in figure [3j In figure [3] we 
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see that the relation in log coordinates suddenly changes slope at T e jf ~ 4000K, the same effective 
temperature where the slope of T e ff-m cz relation (over-plotted in dashed line for ready reference) 
changes. This ma rks the convective cou p ling b etween convection zone and degenerate core, which 
is documented in iFontaine fc van Horn! (|1976l ). Superimposed on the same figure is the relation 
between T e ff and the temperature at the base of convection zone, T b . It is clear that once the 
convection zone couples with the degenerate core, the temperature difference between T c and T b 
becomes approximately constant. 

To understand this phenomenon, let us note the opacity of stellar material can be written in 
the form 

K = K P n T- n - s (1) 

where kq, n and s are in general variables. However, if they approach constant at the base of 
convection zone (which is the case whe n the convection zone couples to the degenerate core-opacity 
in the core is mainly conductive, c.f., IFontaine et al.l ()1974l )). it can be shown that the difference 
in log T between the base of convection zone (represented by the subscript b) and the stellar core 
(represented by the subscript c) is only a function of the temperature gradient at the base of 
convection zone (which equals to the adiabatic gradient V af j as required by the Schwarzschild 
criterion) as follows: 
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In the second line we have changed the upper limit of integration from logP c to 00 because the 
integrand is already vanishing small at P ~ P c (manifested by the fact that log T — log P relation is 
horizontal as shown in figure [1]). The whole integration is done in the conductive/radiative region 
so the formula of radiative temperature gradient is used. We represent the radiative/conductive 
temperature gradient (i.e., ^}°gp ) by V and Taylor expand logT at the base of convection zone: 

1 dV 



log T = log T b + (log P - log P b ) V 



log P=log P b 



p . + (log P-log P b 



2d\ogP 



+ ... (6) 



log P= 



We also note that all derivatives of V with respect to pressure (as well as that of any functions of 
V) at the base of convection zone are functions of V a d| log p =log p fc . (For compactness, we hereafter 
denote the adiabatic gradient evaluated at the base of convection zone simply by V a d, dropping 
the evaluation notation whenever the context is clear): 
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dlogP 
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I log P=log P b 
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(8) 
(9) 



Switching the integration variable to 



e=io g p-iogp 6 



It follows that 



POO 

AlogT« / V ad 10^> v ^d£ = G(V ad ) 

JO 



(10) 



(11) 



This explains the fact that A log T is only a function of V aa \ Since the V aa - at the base of convection 
zone is approximately constant along the cooling sequence, so is AlogT. Thus, in the convective 
coupling regime, we are able to relate T e tt and T c through the following formula: 



logT c -logT ? 



ps 



log P ps 



V ad d\ogP + G(V ad ) 



;P=lo R P h 



(\ogP b -logP ps ){V ad ) + G(V 



ad ) 



(12) 
(13) 



Where the subscript ps stands for the (base of) photosph ere. T px and P ps as functions of T e // are 
found numerically through the radiative transfer code of iHansenl (|1999l ) and is shown in figured 
All of our models show high efficiency of convective transport in the envelope and one need not to 
distinguish between the convective temperature gradient with V ' ad . 

Figure [5] plots T e ff-T c relation for both hydrogen and helium model on the same figure. We see 
that at T e ff > 12 000.fir, the two r elations are almost identical , demonstrating the "radiative-zero 
convergence" (c.f. iKoesterl (119761 ): iFontaine fc van Horn! (119761 )) when the base of convection zone 
and the conductive core is buffered by a radiative zone for both models. In fact, we have already 
seen an indication of this phenomenon in the upper panel of figure [TJ The T e ff—T c relation is thus 
divided into two different regimes: A high T e // regime where T c is insensitive to the atmospheric 
chemical composition, and a low T e ff regime where equation [12] starts to hold for models with 
core-coupled convection zone. 

We are now equipped with the knowledge of atmosphere of two extreme surface composition. 
In the following sections, we will address the complications due to convective mixing. 



4. Self-consistent Treatment of Convective Mixing Envelopes 

As we have mentioned in § [2j most calculations in the current literature assumes the hydrogen 
in post-mixing white dwarf envelope to be highly diluted and has negligible influence on envelope 
structure. As such, both the opacity and equation of state are taken to be that of helium. This 
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assumption is not self-consistent in the sense that it forces the total hydrogen mass (mjj) to go 
from a finite value to zero after convective mixing. Considering that the opacity of hydrogen is 
much larger than that of helium around the T e ff of interest, this assumption needs to be tested. 
Here we take a more rigorous approach to ensure self-consistency in hydrogen content. 

Instead of assuming the chemical composition to go from pure hydrogen to pure helium, we 
assign the post-mixing white dwarf envelopes to possess a (yet unknown) surface composition 
X = X sur f, where X is the hydrogen fraction. 

We take the convective mixing white dwarf atmosphere to be a two layer structure, i.e., a 
homogeneous layer of convective H-He mixture (with X = X sur f) on top of a helium layer, 
which is a direct generalization of the DA white dwarf structure under the influence of convective 
instability (c.f. § 1). The rationale lies in the fact that convective mixing is by far the strongest 
counter-separation mechanism, and diffusion cannot possibly separate the stellar material within 
its range of influence. Besides, we only consider white dwarfs with 10 _11 Mq < mu < 1O~ 7 M0 
to avoid the possible com plications caused by clas sical diffusion (such as prevention and/or delay 



in convective mixing, c.f. (jVauclair &: Reissdll977l )). Thus, the chemical structure of a convective 



mixing white dwarf is described by a step function: 

x=x -"" ,fP< £f7sfe <14) 

where the location of the step (i.e., chemical composition inhomogeneity interface) ought to be 
matched with the boundary of convection zone. 

To proceed, we calculate the envelope structure on a grid of T e ff and X sur f for a given m# 
with the routine described in § [3l which features detailed treatment on the boundary condition 
and EOS of H-He mixture. The resultant models (figure [6]) have structure similar to that of 
non- mixing envelopes (i.e., possess surface convection zone on top of radiative/conductive zone), 
however, the boundary of surface convection zone generally mismatch with the chemical compo- 
sition interface (which is designated by equation [T5] prior to the envelope integration). Our task 
of searching candidates of convective mixing envelopes is thus tantamount to picking envelope 
parameters (m#, Mwd, X sur f, T e *f) which allows a matching between the interface of chemical 
composition and the boundary of surface convection zone. 

A set of concrete examples is provided in figure In line (a) of figure Owe solved the structure 
of a O.6M white dwarf with m fl = 4.0 x 1(T 8 M , X surf = 0.005 at T eff = 4500K. The line 
in the plot depicts the pressure-temperature relation in that white dwarf envelope. The dotted 
line denotes the convection zone while the solid line denotes the radiative zone. According to our 
prescribed X sur j and m#, the interface between surface layer and underlying helium layer locates 
at logP = log x H j inRi = 17-3) which is indicated by the long vertical dashed line. We found 
that this model cannot be a candidate of convective mixing envelope because while the layer of 
X sur f = 0.005 extends to logP = 17.3, the convection zone merely reaches logP = 17.18. This 
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configuration cannot be realized because hydrogen in such envelope would diffuse upwards, altering 
the prescribed chemical profile. In other words, the convection zone mass in the resultant model 
underestimates the self-consistent convection zone and we do not consider this model as a candidate 
of self-consistent, convective mixing envelope. 

In line (b) of figure [6] we solved white dwarf atmosphere with identical run and X sur f, but 
T e fj = 3800-KT. The resultant model cannot be realized as well because the convection zone now 
extends to logP = 17.44, which is below the H-He mixture/ He layer interface. The convective 
motion would mix pure helium gas with the gas of X = 0.005, resulting in a gas which is more 
helium-rich. In this case, the convection zone mass in our model overestimates that of self-consistent 
models. However, a self-consistent model can be found between the two mentioned cases with 
T eff ~ 4128K (line (c)). 

In mathematical terms, the parameter that allows self-consistent convective mixing envelopes 
is defined by the contour: 



We have calculated the contour for a wide range of run and Myy £>. Some representative 
contours are plotted in figure [7] (For the rest of this paper, such contour is referred as a "contour 
of hydrogen conservation"). We clearly see that the possibility for convective mixing envelope to 
possess a non-negligible surface hydrogen fraction. Therefore, the assumption of highly diluted 
surface composition in convective mixing envelopes is not always proper. 

To conclude, we have placed a constraint on the possible outcomes of convective mixing: 
Namely, the post-mixing envelope must lie on the contour of hydrogen conservation in order to 
conserve the total hydrogen mass. In the following sections, we will exploit this result and derive 
the evolutionary sequence, cooling curves and chemical evolutionary curves for convective mixing 
white dwarfs. 



We now formally revisit the convective mixing scenario in white dwarfs with intermediate 
amount of hydrogen. As we have reviewed in § [21 the surface convection zone in a hot hydrogen 
white dwarf deepens continuously as T e ff decreases. Its spectrum remains DA until the lower 
boundary of convection zone reaches the H /He boundary. At this point, convective mixing occurs 
and the X sur f is bound to change. We claim the outcome is governed by two rules: 

1. The total amount of hydrogen must be conserved. This rule can be met by requiring the 
post-mixing envelope to lie on the contour of hydrogen conservation. 

2. T c is approximately invariant. The underlying principle of this rule is energy conservation. 
White dwarf is a highly degenerate system, therefore, the gravitational potential energy re- 




surf 



(16) 



5. Revisiting the Convective Mixing Scenario 
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Fig. 1. — Pressure-temperature relation of four exemplary white dwarf envelope models. Convection 
zone is represented by dotted lines and radiative/conductive zones are represented by solid lines. 
The models with hydrogen surface are marked by "(a)" while those with pure helium surface are 
marked by "(b)". The vertical dashed line indicates the interface between H and He in hydrogen 
surface models and has no meaning pertaining to the helium models. Upper: The P-T relation at 
logT e ff = 4.3. Lower: The P-T relation at logT e ff = 3.8. It is clear that the core temperature 
T c is relatively insensitive to the surface boundary condition at high T e ff, manifesting so-called 
"radiative-zero convergence" . 



10 4 




3000 6000 9000 12000 15000 



Fig. 2. — T e ff-m cz relation for a O.6M white dwarf hydrogen envelope. The convection zone 
deepens with the decrement of T e ff until T e ff ~ 4000i"C. The maximum of convection zone mass 
is around 1.5 x 10~ 6 Mq for a O.6M0 white dwarf. 
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Fig. 3. — T e ff-T c and T e ff-Tf, relation for a 0.6M Q hydrogen white dwarf envelope. Over-plotted 
on the same figure is T e ff-m cz relation (dashed line), which we have seen in figure EJ We can see 
that the T e ff where the slope of T e jj~T c relation changes coincides with the T e ff where the slope 
of T e ff-m cz relation flattens (both slope measured in log coordinates). This T e ff is marked by the 
long vertical dotted line in the figure. 
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Fig. 4. — The T e ff-P ps relat ion and T P ff—T ps relation of pure hydrogen photosphere obtained from 
the radiative transfer code of lHansenl (|1999l ). Upper: T e ff 
The dotted line is the Eddington T-r relation T 
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Fig. 5. — T e ff-T c relation for DA(solid) and DB (dash-dotted) white dwarfs. 
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Fig. 6. — Three white dwarf envelope models with X sur f = 0.005 and ran = 4.0 x 10 M . 
The chemical composition interface is determined by X sur f and ran and is denoted by the long 
vertical dashed line. The three models have different T e ff. In (a), T e ff = 4500i"C. The resultant 
convection zone does not reach the interface, yielding an inconsistency. In (b), T e ff = 3800-ftT. The 
resultant convection zone penetrates the interface, yielding another inconsistency. However, in (c), 
T e ff = 4128-ftT, the boundary of surface convection zone and the chemical composition interface 
match. The conclusion we draw here is the following: If there is a white dwarf with X sur f = 0.005 
and m H = 4.0 x 1O~ 8 M , then T e ff can only be 4128K at this neighborhood. 



-12- 



leased during convective mixing is small compared to the total thermal energy, which is 
represented by T c . 

It also follows from the reasoning of principle [2] that T c is monotonically decreasing over time 
(because heat energy is continuously released through radiation). It can therefore serve as a sorting 
index and we proceed to calculate the T c for all models that are possible to realize during the 
evolutionary sequence. This include all pure DA envelopes with T e fj > T m i x (the convective 
mixing T e ff) and the convective mixing envelopes that lie on the contour of hydrogen conservation. 

An exemplary set of results is summarized in figure El where the relations between T c , T e ff 
and X sur f is plotted for ~ 200 candidate models. 

It is important to note, however, that the contour of hydrogen conservation is merely a con- 
straint on possible outcomes. Models on the contours need not to be realized in the cooling sequence. 
Let us denote the core temperature of the coolest pure DA star (among the candidate models) to 
be T*. This particular model, as we have reviewed, is just about to be convectively-mixed. We 
also know that prior to this model, convective mixing is not present and the white dwarf has pure 
hydrogen surface, due to the efficient separation of diffusion. Therefore, all models on the contour 
of hydrogen conservation (i.e., those does not have pure hydrogen surface) with T c > T* should be 
excluded from the cooling sequence. 

The evolution of white dwarf in T e f f and X sur f is thus found by sorting the models according 
to T c after the exclusion of unrealized models. It is indicated by the arrows in upper and middle 
panels of figure where we have marked the excluded models in different symbol. In our example 
(tuh = 10~ 8 M Q ), it is apparent that X sur f decreases and T e ff increases upon convective mixing. 
The evolution track of white dwarf in T e ff-X sur f space is therefore zigzag, as shown in the lower 
panel of figure [8l The white dwarf travels through a segment of T e t j interval twice during its 
course of evolution. The first time as a pure DA star, the second time as a white dwarf with lower 
X SU rf- Incidentally, we have examined the behavior of mH(T e ff,X sur f) to ensure the stability of 
our solutions. 



6. The Cooling curves and the Chemical Evolution curves 



In order to obtain the cooling curve (t-T e f f relation) and the chemical evolution curve {t-X sur f 
relation), one needs not only to know the evolution track in the parameter space (that constructed 
by T e ff, T c and X sur f) but also needs to know the time spent between different evolutionary stages. 
This sets the goal of this section. 



With the help of a full evolutionary model of white dwarfs ([Hansen! 119991 ) , we are able to 
calculate the energy of a white dwarf a s a function of its core temperature, i.e., E = E(T C ). 
The evolutionary model of IHansenl (119991 ) is designed to deal with models of a fixed atmospheric 
composition. However, it is sufficient for our purposes since the energy content in the atmosphere 
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is negligibly small. Our result is shown in figure [H 

The cooling time of a particular evolutionary sequence can be calculated by relating T e ff to 
the core temperature and energy content: 



dt 



dE(T c ) 
' 4xR*oT*.JT c ) 



t(T c 



E{T C ) 



dE 



4 * R2 ° T eff 



dE(T c ) 
dT c 



dT, 



4*R*<jT? ff 



(17) 
(18) 



Note that E(T C ) is scenario independent whereas T e ff(T c ) depends on X sur f and is specific to 
the particular evolutionary sequence. 

We are now able to calculate the chemical evolution curve and T e ff evolutionary curve in the 
following way: 



1 . For a given run , we first identify the evolutionary sequence of a white dwarf in the convective 
mixing scenario. T e ff(T c ) and X sur f(T c ) is thus obtained. This process is described in 
section [5] and an exemplary result was shown in figure 

2. A fter the e voluti onary sequence is identified, we use E(T C ) from the full evolutionary model 
of iHansenl (|1999l ) and T e ff(T c ) to calculate the cooling time t(T c ) through equation [TBI 

3. T e ff(t) and X sur f(t) are then obtained by inverting t(T c ) and substituting it into T e ff(T c ) 
and X surf (T c ). 



A few examples of these curves are shown from figure [TU] to figure [TH 



7. Discussion 



In this section, we discuss the physical meanings of the cooling curves. We also relate our 
results to the reported observations. 



7.1. Comparison Between Models 

We would like to compare the cooling curves for pure hydrogen stars {mjj = 1O~ 5 M0), pure 
helium stars, and stars that undergo spectral evolution (here we use a 0.6Mq white dwarf with 
nriH = 1CP 8 M as a representative). Strictly speaking, there are two types of cooling curves. In the 
literature one usually refer "cooling curve" to the t—T e ff curve (as we did in the previous sections). 
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In this section, however, we will reserve the term "cooling curve" for the core cooling curve, t—T c 
relation. t—T e fj relation will be temporarily termed "fading curve". This difference in terminology 
is made for the sake of clarity, as we will see below. 

A comparison of the fading curves are shown in the upper panel of figure [15l The fading curve 
of chemical-evolving white dwarf (mu = 10~ 8 Mq) is seemingly complicated, however, in essence 
it is the relatively simple cooling curve (as shown in the lower panel of figure fT5|) masked by the 
(discontinuous) T e ft—T c relation which is given by figure [T6j The cooling curve in the lower panel 
of figure [15] shows that the cooling rate is generally higher for the white dwarfs with lower X sur f 
(fixing T c ). 

The cooling rate as a function of T c is directly determined by the relation between T e ff and 

T ■ 



v ; dt dt dT c y 1 

^ _ _ ^ R2T eff(T c ) 

J4- ~ dE ~ dE \ ZU > 

Thus, models with closer T e tf—T c relation have higher cooling rates (fixing T c ). T e ff—T c relation for 
all three models in the convective coupling regime is shown in figure [TBJ The triangular symbols are 
the relation for chemical-evolving models and the two bracketing curves are the relation for pure 
hydrogen stars (upper) and pure helium stars (lower). We can see that in the regime of our interest, 
T e ff~T c relation is significantly different for pure hydrogen models and pure helium models. That of 
chemical evolving white dwarfs, on the other hand, varies between the two extremes. Of particular 
interest is the discontinuity of T e ff-T c relation located at logT c = 6.53 for chemical evolving white 
dwarfs, which corresponds to the onset of convective mixing. The mixing leads to a sudden increase 
in helium content and T e ff. 

The mechanism of increase in T e ff is a combination of convective coupling and reduction of 
surface opacity. Note that both ingredients are important. In the complete absence of convective 
coupling, T c is insensitive to photosphere opacity. It is only after convective coupling, a reduction 
in surface opacity allows corresponding reduction of log T c -log T e f t through the increment of P ps 
in equation [T2l 

Fading rate, on the other hand, is determined by both cooling rate and the derivative of T e ff 
to T c : 

dT eff = dT c dT eff 
dt dt dT c K ' 

As a result, a closer relation between T e ff and T c does not imply higher fading rate. For example, 
the fading rate of helium models at t < lGyr in figure [15] is lower than that of hydrogen models, 
despite of their higher cooling rate. 
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7.2. Comparison with Previous Works 



This work is the first study on convective mixing that accounts the post-mixing composition in 
a self-consistent fashion. B esides this , howe v er, there are additional di ffere nces between th is study 
and previous works such as IShipmanl (|1972l ) , iBaglin k, Vauclairl (|1973l ) and iKoesterl ()1976l ) . 



Shipmanl ()1972l ) based his argument of T e f f increment solely on the reduction of opacity, which 
is not complete. Convective coupling is crucial for T e ff to increase upon convective mixing. In 
particular, the increment of T e ff in the example that Shipman used (13000-KT DA star) turns out 
to be negligible in our calculation due to the lack of convective coupling. 



Baglin &: Vauclairl (|1973l ) assumed the physical conditions below the post-mixing convection 
zone to be invariant during the mixing process, which is not accurate. Our analysis shows a 
non-negligible change in logT at the mass fraction which corresponds to the base of post-mixing 
convection zone (c.f. figure fT7|) . We note that the luminosity is constant throughout the envelope 
and fixing the physical conditions at any location other than the isothermal core will naturally fix 
the total luminosi ty. Beside s , the authors claimed that the inclusion of convection zone reduces the 
T e ff increment of IShipmanl (|1972l ). contrary to our findings. 



The result of lKoesterl (jl976l ) (no luminosity change upon convective mixing) is different from us 
because he only investigated one particular model with extremely small mn- Due to the smallness 
of thh, the convective mixing occurs at very high T e ff where T e ff—T c relation is still insensitive to 
surface composition (see section [3]). 



7.3. Spectral Evolution 



We have successfully calculated the evolution of X sur f for white dwarfs with m# varying over 
five orders of magnitude. Howeve r, the evolution of spectra is not quite the same as the evolution 



of X sur f. 
~ 12000i<C 



Bergeron et al.1 (I199Q ) stated that helium lines start to turn invisible at T e ff below 
As a result, white dwarf atmosphere with log ^y^gy < 2 would all appear to be of 



spectral type DA, even if in reality, helium could be the main constituent. 



In lBergeron et al.l (j2007l ). it was again emphasized that the spectroscopic mass of low T e ff DA 
white dwarfs has a much greater mean and dispersion than its high T e ff counterparts. The authors 
thus speculated that upon convective mixing (i.e., when hydrogen convection zone breaks into the 
helium layer), white d warf stars might still appear as spectral ty pe D A instead of turning in t o tha t 
of spectral type DC as IShipmanl (|1972l ). lBaglin fe Vauclairl (|1973l ) and lFontaine fc Wesemaell (|199ll ) 
have proposed. 

Our results show that both proposal has its own standing. In our higher mu models (mjj ~ 
10~ 7 M Q ), the white dwarfs adjust their X sur f to O(0. 1) and T e ff increases only slightly (~ 
O(100K)) upon convective mixing. They will thus appear to astronomers as a slightly hotter 
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and spectroscopically more massive DA white dwarf. 

On the other hand, DA white dwarfs with ma < 10~ 9 Mq will directly turn into DC stars 
upon convective mixing. The T e // wm increase at the order of lOOOi'T. 

What we have shown is, upon convective mixing, a white dwarf star always decreases its X sur f. 
However, the spectral outcome and T e tf increment will be dependent on run- 



8. Conclusions 

In this paper, we have developed a scheme to calculate the possible atmosphere configurations 
of a white dwarf, given its stellar mass (Mvkd) and hydrogen content Based on the knowledge 

of these possible configurations, we worked out the cooling curves and chemical evolution curves of 
a white dwarf of fixed mj, undergoing convective mixing. 

We would like to emphasize that although the main results are only given in a convective 
mixing scenario in the absence of accretion, our framework is actually versatile and is not confined 
to such scenario. The result in section [J] encompasses hydrogen content over multiple magnitudes 
and is applicable to the case where hydrogen content is variable (either due to wind or accretion). 
To calculate the cooling curve in a scenario where mn is variable, one could pick the stages from 
the larger set of possible configurations (with many different m#). Once the configurations are 
picked one could calculate the cooling curve and chemical evolution curve in similar fashions. 

We have adopted a semi-analytic approach to comprehensively study the possible configura- 
tions in the pr oblem. We have confirmed our results with full white dwarf evolutionary code (e.g., 



(jHansenl Il999l )) by switching X sur f from 1 to the value we obtained from this work. We found 



th at T P ff will au tomatically increase accompanying the decrease in X sur f, contrary to the results 



of iKoesterl (119761 ). The flexibility of our semi-analytic models enables us to examine the full range 
of behaviors associated with different hydrogen layer masses, thereby encompassing a variety of 
scenarios presented in the literature. 

In a later paper, we will apply these cooling curves and X sur f evolutionary curves to calculate 
non-DA to DA ratio as a function of T e ff, in attempt to solve the problem of the "non-DA gap" 



(jBergeron et al.lll997l ). We will also obtain luminosity functions from them and discuss the impact 



of chemical evolution on cosmo-chronology. 



The work described here is supported by NASA grant ATP03-000-0084 and the Alfred P. 
Sloan Foundation. E.C. would like to thank Dr. Bernard Freytag for discussions on the subject of 
convective overshooting. 
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Fig. 7. — Upper: The contour of hydrogen conservation for a 0.6M© white dwarf with various ran 
(value marked on the contour line). Lower: The contour of hydrogen conservation for white dwarf 
of various mass (value marked on the contour line) with mn = 10 _9 M©. 
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Fig. 8. — T c -T e ff relation (upper panel), T c -X sur f relation (middle panel) and T e jj-T c relation 
(lower panel) for white dwarf envelopes with a total hydrogen amount of 1O~ 8 M from self- 
consistent, convective mixing calculation. The models which will be realized during the cooling 
sequence are marked by triangular symbols while those will not are marked by hollow circles. The 
arrows indicate the track of evolution, based on T c . We can see in the lower panel that the track 
is zig-zag in T e ff-X sur f space. 
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log T c 



Fig. 9. — The relation between logT c and thermal energy. The zero point of energy is set at the 
beginning of the evolution sequence, where logT c = 7.7 (hence the negative value of energy). 




t(Gyr) t(Gyr) 

Fig. 10. — tcooiing y s T e ff and t coo u ng vs X sur f for white dwarf envelopes with run = 1-0 x 10 M Q . 
The time spent between each white dwarf configuration is calculated by equation 1181 We can see 
here that a white dwarf of mjj = 1.0 x 1O~ 7 M0 will probably remain as a DA after dredge-up, 
although its X sur j actually decreases by about two orders of magnitude. The T e ff does not decrease 
monotonically with time. Upon dredge-up, its value increases (~ O(lOO-fT)). 
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Fig. 11. — Cooling curve and Chemical Evolution curve for white dwarf envelopes with ran = 
1.0 x 10~ 8 M Q . 
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Fig. 12. — Cooling curve and Chemical Evolution curve for white dwarf envelopes with ran 
1.0 x 1O~ 9 M . 
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Fig. 14. — Cooling curve and Chemical Evolution curve for white dwarf envelopes with mn = 



1.0 x 10 _11 M f 
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Fig. 15. — The fading curve (upper panel) and cooling curve (lower panel) of 0.6Mq white dwarf 
with three different m#: uih = 1-0 x 1O~ 5 M0 (dashed line), mjj = 1.0 x 10~ 8 Mq (solid line) 
and mn = (pure helium white dwarf, dash-dotted line). The relation between cooling rate and 
hydrogen content is fairly simple: white dwarf with lower X sur f cools faster. However, the relation 
between fading rate and hydrogen content is more complicated and is given by equation [2TJ 
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Fig. 16. — T e ff—T c relation for white dwarf models in figure [T5l That of the thick hydrogen model 
is indicated by the dotted line and that of pure helium model is indicated by the dashed line. 
The relation of the convective mixing models is represented by the triangular symbols, which lies 
between the two non- mixing cases. 
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Fig. 17. — The pre-mixing model (X sur f = 1, T e ff ~ 7967i"T, labeled by "(a)") and post-mixing 
model (X surf ~ 0.0003, T eff ~ 8750i^T, labeled by "(b)") of a 0.6M© white dwarf with m H = 
1O~ 9 M . We have shown that the temperature at the mass fraction which corresponds to the base 
of post-mixing convection zone is changed during the convective mixing process. 



